Two-gap superconductivity in single crystal Lu2Fe3Si5 from penetration depth 

measurement s 
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Single crystal of Lu2Fe3Si5 was studied with the tunnel-diode resonator technique in Meissner 
and mixed states. Temperature dependence of the superfluid density provides strong evidence 
for the two-gap superconductivity with almost equal contributions from each gap of magnitudes 
Ai/fcsTc = 1.86 and Ai/fcsTc = 0.54. In the vortex state, pinning strength shows unusually 
strong temperature dependence and is non-monotonic with the magnetic field (peak effect). The 
irreversibility line is sharply defined and is quite distant from the Hc'ziT), which hints on to enhanced 
vortex fluctuations in this two-gap system. Altogether our finding provide strong electromagnetic 
- measurements support to the two-gap superconductivity in Lu2Fe3Si5 previously suggested from 
speciflc heat measurements. 

PACS numbers: 74.25.Nf, 74.25. Op,74.20.Rp,74.25.Ha 



Originally, the interest to iron - containing silicides 
M2Fe3Si5 was due to unusually high superconducting 
critical temperatures for compounds containing a crys- 
tallographically ordered iron sublattice 1]. Three sili- 
cides, all having the same tetragonal structure, have been 
found to be superconductors, M=Y, Sc and Lu with tran- 
sition temperatures of 2.4, 4.5, and 6.0 K, respectively. 
It turns out that iron in these materials is nonmagnetic 
as was concluded from ^^Fe Mossbauer effect measure- 
ments 0, However, further detailed studies revealed 
that other superconducting properties are quite uncon- 
ventional. The upper critical field Hc2(0) for Lu2Fe3Si5 
has been found to be unusually large when compared 
to other iron-containing superconductors 0, Q and its 
temperature dependence differs from conventional. The 
anisotropy and a pronounced peak effect in magnetic 
measurements was reported in 17]. The presence of a 
large residual electronic term in the specific heat below 
Tc as well as a reduced specific heat jump at Tc have been 
observed and confirmed, indicating departures from the 
standard BCS-like behavior [^, Q . Non-magnetic impuri- 
ties suppressed Tc at a significant rate, incompatible with 
isotropic s-wave BCS picture [1,0]. On the other hand, 
ac Josephson effect indicated s-wave pairing mechanism 
[lo| . Vining et al. have proposed a two-band model 
in order to explain their specific heat data Q. Their 
model assumes two-band Fermi surface with one band 
being superconducting and gapped, and another being 
normal. This represents an extreme case of multiband 
superconductivity as we know it today, for example in 
MgB2 superconductor where different bands have gaps 
of different magnitude 11 , 12] • Later detailed measure- 



ments of Lu2Fe3Si5 crystals and analysis have shown that 
specific heat data are well explained quantitatively within 
two band model of superconductivity where both bands 
are gapped but with different gap amplitudes [13] • 

In this letter we present precision measurements of 



the London and Campbell penetration depths, analyze 
superfluid density as well as unusual vortex properties 
and conclude that Lu2Fe3Si5 is, indeed, a two-gap super- 
conductor. It seems that multiband superconductivity is 
more widespread and develops when there is different di- 
mensionality of the Fermi surface on different bands that 
leads to a reduced interband scattering. In MgB2 there 
are two- and three - dimensional bands [ll|, [l3] , whereas 
Lu2Fe3Si5 has quasi-one- and three- dimensional Fermi 
surfaces |13| . 

Measurements of Lu2Fe3Si5 single crystal were per- 
formed using a tunnel diode resonator (TDR) [3, iM l3]- 
Extended review of using TDR to study superconductors 
is given in Ref.[l3]. The main components of the TDR 
are an LC tank circuit and a tunnel diode. The tunnel 
diode has a region of negative differential resistance in its 
I-V curve. If a DC bias voltage is applied across the diode 
in this region, then it acts as an AC power source for the 
LC tank circuit. This results in a self-oscillating circuit, 
which resonates continuously at a constant frequency for 
given values of L and C. The resonance frequency of the 
circuit used in our measurements was near 14 MHz. All 
throughout the measurements the circuit is kept at a con- 
stant temperature, 4.8 K ± 1 mK, allowing for a stability 
of 0.05 Hz in the resonance frequency over several hours. 
The sample to be studied is mounted on a sapphire rod 
with a small amount of Apiezon N grease. The sapphire is 
inserted inside of the inductor coil of the tank circuit. It 
is important that the sample and its mount do not make 
physical contact with the coil so that the temperature of 
the sample may be changed while keeping the circuit at 
a constant temperature to maintain the stability. As the 
magnetic susceptibility of the sample changes with tem- 
perature, so does the inductance of the tank coil. This 
results in a change in the TDR resonance frequency. By 
measuring the shift in the resonance frequency, we are 
able to sense changes in the penetration depth on the 
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order of 0.5 Angstroms. Specifically, the frequency shift, 
A/ = / (T) — /o, with respect to the resonant frequency 
of an empty coil, /o, is given by 



A/ (T) = -GAttx (T) = G 



A , fR 



(1) 



where G ~ /oFs/2Vc (1 — N) is the geometry dependent 
calibration constant, Vs is sample volume, Vc is the ef- 
fective coil volume and N is the demagnetization factor. 
The effective sample dimension R is calculated by using 
Ref . [isj . As described in detail in Ref.[l^ it is difficult 
to obtain the absolute value of the penetration depth 
due to uncertainties in the sample dimension. How- 
ever, it is possible to calibrate the system with great 
accuracy by using temperature-dependent skin depth, 
S (T), measured right above the Tc. In that regime, both 
real and imaginary parts of the susceptibility are taken 
into account and the frequency shift due to skin effect 
is A/(r)y>y^ = G[l - ((5/2i?)tanh(2i?/(5)]. The skin 
depth, 5 (T) = cy^ p (T) /2Truj is evaluated independently 
from the temperature-dependent resistivity, p{T), mea- 
sured by the four-probe technique. In addition to excel- 
lent stability and sensitivity, the advantage of this tech- 
nique is very low excitation fields, ^ 20 mOe which en- 
sures that the sample is in Meissner state. Furthermore, 
by superimposing an external DC field we can probe the 
vortex state in so-called Campbell regime where small 
excitation ensures that vortices remain in their potential 
wells. 

Single crystal of Lu2Fe3Si5 was grown by the floating- 
zone technique using an image furnace followed by an 
annealing as described in detail elsewhere [isj . The sam- 
ple was a rectangular slab having dimensions 0.99 x 
0.84 X 0.15 mm'^ with the c-axis perpendicular to the 
largest face. To study possible anisotropy of the response, 
the measurements were performed for the excitation field 
both parallel and perpendicular to the c-axis of the sam- 
ple. A "^He cryostat with sample in vacuum and external 
field up to 9 T was used for the reported studies. 

Figure [1] shows the temperature dependence of the 
London penetration depth, Xab(T) obtained from the 
measurements along and perpendicular to the c-axis. 
Both orientations give XabiT), because sample is a thin 
plate and apparently Ac(T') is not too different from the 
Xab{T) - otherwise the results would not coincide. The 
value of Aab(O) = 0.2 pm was obtained as described 
in Ref. [31 from the reversible magnetization dM/dlnH 
measured independently on the same sample using Quan- 
tum Design magnetometer. In the further analysis, pos- 
sible uncertainty of this number up to 25% was examined 
and and confirmed not to change our conclusions in any 
way. 

Symbols in Fig. [2] show temperature dependent super- 
fluid density, ps{T) = (A (0) /A {T)f calculated from the 
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FIG. 1: (Color Online) XabiT) obtained from the measure- 
ments in two orientations. Inset: Evaluation of the calibra- 
tion constant by matching Atix{T) from the TDR data (empty 
points) to calculated from the the skin depth (solid red curve) 
obtained from the resistivity. 



penetration depth. Fig. [TJ The solid red curve is the to- 
tal superfluid density calculated from the a model that 
assumes two independent contributions to the total su- 
perfluid density and has been successfully applied to the 
well known two-gap superconductor MgB2 [ll|, . In 
this model, each superconducting gap, Ai(T) and A2(T) 
have similar temperature dependence given by the weak- 
coupling BCS self-consistency equation j but with dif- 
ferent ratio of /S.{Q)/kBTc that become two fit parame- 
ters. A third fitting parameter gives the relative con- 
tribution of each band to the total superfluid density, 
Ptotai{T) = xpi(T) + (I-x)p2(r). Each superfluid den- 
sity is calculated by using full temperature range semi- 
classical BCS treatment as described in detail elsewhere 
14| . These partial pi (T) and p2 {T) are shown by marked 
solid lines in Fig. [H The best fit was achieved with 
X = 0.51, Ai/zfesTc = 1.86, and Ai/ksTc = 0.54. The 
first gap is quite close to the weak-coupling value of 1.76, 
whereas the second gap is much smaller and it is surpris- 
ing that earlier two-band model assumed it to be fully 
normal [6|]. Similarly to MgB2, the two gaps contribute 
equally to the superfluid density. A dashed line, which 
almost follows the data is calculated from the param- 
eters obtained analyzing specific heat data, x — 0.47, 
Ai/keTc = 2.2, and Ai/ksTc = 0.55 [Hi, which is in 
a quite good agreement given very different nature of 
the measurements. To further highlight the qualitative 
difference between single and two gap behavior, we plot 
dps/dt in the inset to Fig. [21 Note characteristic non- 
monotonic behavior in the case of two gaps. It is not 
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FIG. 2: (Color Online) TDR data (symbols) fitted to a two- 
gap model with indicated parameters. The red curve shows 
the total superfluid density. Also shown and labeled two par- 
tial superfluid densities as described in the text. The dashed 
line is calculated fro the parameters derived from the specific 
heat data. Inset: dps/dt for pure d-wave, s-wave and present 
case of two-band superconductivity. 



present either in pure d-wane, nor pure s-wave cases. 

Whilst the situation is quite clear for the London pene- 
tration depth, measurements in an applied magnetic field 
reveal more puzzling behavior of the studied compound. 
When an external DC field is applied and small - am- 
plitude AC response is probed, vortices respond elasti- 
cally and the overall susceptibility is governed by the 
Campbell penetration depth, = A| + A^, where Ai, 
is the usual London penetration depth described above 
and Xc{B,T, j) is the Campbell penetration depth 

^ 4>viB/A'Ka{i). Here 0o is the flux quanta and a{j) 
is the Labusch parameter that generally depends on the 
biasing Bean current generated in the sample, for exam- 
ple, after applying field after cooling in zero field. Sample 
magnetic susceptibility (and the frequency shift) in the 
vortex state is still given by Eq. [l] but with generalized 
penetration depth. 

In conventional type-II superconductors there is no 
hysteresis for zero-field cooled (zfc) and field cooled 
(fc) curves of the small amplitude AC response. How- 
ever, in materials where jc is strongly temperature de- 
pendenl] (e.g. high- Tc cuprates) large hysteresis is ob- 
served 



17|. As shown in Ref.[l3|, cubic correction to 
a parabolic potential well for vortex pinning leads to 
=ao\/l -j/jc, where jc = cagrp/^o is the criti- 
cal current and Vp is the radius of the pinning potential. 
This model explains why zero-field cooled curve differs 
from subsequent cooling and warming and it was suc- 
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FIG. 3: (Color Online) 47rx(r) from the TDR data at three 
diflerent values of an applied magnetic field along the c-axis. 
Each curve was obtained after cooling in zero field and then 
warming and cooling twice. The labels and arrows indicate 
various characteristic points used in later analysis. 



cessfuUy used to explain the data for Bi2Sr2CaCu208+y 
superconductor. 

47rx(T') in the vortex state of Lu2Fe3Si5 is shown in 
Fig. [3] for three representative fields. In each case, sam- 
ple was cooled in zero applied field to the base temper- 
ature and indicated magnetic field was applied. Then 
measurements were taken while warming up the sample 
above (zfc-w). Then, sample was cooled and warmed 
twice without changing the field and while taking the 
data (fc-c and fc-w). For low field values, there is no 
hysteresis observed, while at intermediate fields the hys- 
teresis becomes very pronounced. Clearly, the hysteresis 
is associated with the static Bean current, j, induced by 
applying field. We also note that this effect is not associ- 
ated with the vortex density (e.g. less vortices after zfc), 
because then initial Campbell length would be smaller 
than equilibrium, not larger as observed. 

By measuring many 47rx(r) curves at different mag- 
netic fields, we extracted field dependence of the initial 
susceptibility obtained after zfc and fc. Figure H] shows 
the resulting AirxirriH) (open circles) and 4:TTXrev{H) 
(closed squares) curves at T = 0.7 K. The inset shows 
the difference between the two curves. This difference 
is directly related to the strength of pinning and magni- 
tude of the apparent Bean current density, j, Ax ~ j / jc 
where we assumed j <^ jc- There is a clear peak effect 
and its location is well compatible with direct measure- 
ments reported in Ref. 

Finally, we construct the H — T phase diagram ob- 
tained from our measurements for both directions. While 
Meissner response is governed by currents flowing in the 
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FIG. 4: (Color Online) Open circles: 4:Tvxirr{H) at T = 0.7 
K measured by applying field after zfc as indicated in Fig. [S] 
Closed squares indicate 4:TTXrev{H) obtained on fc. Inset: the 
difi'erence between zfc and fc curves shown in the main frame, 

AttAx = 4-KXirr — AlVXrev ■ 




crystal in two orientations. 



ab-plane, in magnetic field the response is anisotropic 
and is determined by orientation of vortices with respect 
to crystal axes. We observe large anisotropy of the up- 
per critical field, Hc2{T), down to 1 K as shown in Fig.[5l 
which has not been reported in earlier papers. Further- 
more, Hc2{T) determined from the TDR measurements is 
in excellent agreement with the specific heat data. Note 
that Hc2{T) is linear in temperature down to O.lSTc. Fig- 
ure also shows position of the irreversibility line (see 
Fig. [3] for definition) for both orientations. Unlike con- 
ventional superconductors where Hirr (T) is very difficult 
to determine, because it gradually merges into Hc2 (T) , in 
Lu2Fe3Si5 it is sharply defined and is quite distant from 
the Hc2{T), which is another indication of significant re- 
duction of the critical current possibly due to enhanced 
fluctuations in the two-gap system. 

In conclusion, we find that Lu2Fe3Si5 shows Meissner 
response compatible with two-gap s-wave superconduc- 
tivity, similar to MgB2. It seems that having Fermi sur- 
faces of different dimensionality is the important ingredi- 
ent for multi-gap superconductivity. In the vortex state, 
Lu2Fe3Si5 shows unusually strong temperature depen- 
dence of the critical current, which is also non-monotonic 
with magnetic field (peak effect). The upper critical field 
is anisotropic and linear in temperature. All these ob- 
servations are reminiscent of unconventional supercon- 
ductivity and further theoretical insight to connect these 
properties is needed. 
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